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一 种 非 线 性 尺度 空间 自 适应 均衡 水 印 算法 
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(辽宁 工程 技术 大 学 软件 学 院 , 辽宁 芯 芦 岛 125105) 


摘 要 : 针对 线性 尺度 空间 水 印 算 法 嵌入 水 印 位 置 定位 不 够 精确 ， 谨 入 强度 参数 选取 随机 ， 提 出 一 种 非 线性 尺度 空间 
自 适应 均衡 水 印 算法 。 利 用 KAZE 算法 提取 并 筛选 出 非 线 性 尺度 空间 稳定 性 强 的 特征 点 ， 构 建 谱 入 水 印 区 域 ; 将 水 印 
图 像 奇 异 值 分 解 ， 构 造 新 给 阵 作为 待 座 入 水 印 载体 ， 通 过 调整 果 蝇 优 化 算法 的 适应 度 函 数 计算 嵌入 强度 ， 结 合 DWT- 
SVD 算法 自 适 应 完成 识 入 水 印 过 程 ,对 受到 攻击 的 水 印 图 像 提 取 特 征 点 , 合成 特征 区 域 矩 阵 , 使 用 瞪 入 水 印 的 逆 过 程 ， 
提取 水 印 。 实 验 结果 表明 ，PSNR 值 均 达 到 44dB 以 上 , 平均 NC 值 高 达 0.99， 有 效 均衡 了 水 印 算 法 的 不 可 见 性 和 鲁 棒 
性 。 
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Adaptive equalization watermarking algorithm in nonlinear scale space 


Qi Xiangming, Li Shuang, Li Yue, Hou Mingjun 
(College of Software Liaoning Technical University, Huludao Liaoning 125105, China) 


Abstract: Aiming at the inaccuracy of embedding watermarking location and random selection of embedding intensity 
parameters, this paper proposed a nonlinear scale space adaptive equalization watermarking algorithm. It used KAZE algorithm 
to extract and filter the feature points with strong stability in nonlinear scale space for constructing the watermark region. Then 
it used singula value decomposition on watermarking image, constructed a new matrix as the watermarking carrier to be 
embedded and calculating the embedding intensity by adjusting the fitness function of the ftuit fly optimization algorithm, and 
with the DWT-SVD algorithm, the watermark embedding process is adaptively completed. Extracting the feature points from 
the under attack watermark image to synthesize the feature region matrix. And it used inverse process to extract watermark 
image. The experimental results show that the PSNR values are above 44dB and the average NC value is as high as 0.99, which 
effectively equalizes the invisibility and robustness of the watermarking algorithm. 
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图 像 变换 ， 再 结合 奇异 值 分 解 ， 提 高 了 嵌入 水 印 的 鲁 棒 性 ， 但 
因 只 和 峙 入 水 印 奇异 值 向 量 ， 没 有 水 印 图 像 基 空 间 结构 信息 ， 不 

大 数据 时 代 每 天 都 会 有 海量 信息 产生 ， 网 络 共享 使 得 数字 可 见 性 较 差 。 文 献 [ 和 将 NSCT、DCT、DWT 和 SVD 相 结 合 一 
产权 保护 显得 越 来 越 重要 。 水 印 作为 数字 产权 保护 的 一 项 有 效 。” 定 程度 均衡 了 不 可 见 性 和 重 棒 性 ， 可 以 有 效 抵抗 常规 攻击 ; 文 


措施 ， 如 何 有 效 解决 图 像 版 权 防 伪 认 证 ， 快 速 批量 处 理 图 像 水 。 献 [5, 6] 主 要 是 为 了 提高 算法 鲁 棒 性 , 但 只 能 抵抗 常规 攻击 , 对 
印信 息 ， 成 为 热门 研究 课题 。 剪 切 、 旋 转 等 几何 攻击 鲁 棒 性 较 差 。 
水 印 是 嵌入 载体 后 不 易 被 觉察 , 且 不 影响 载体 使 用 价值 的 选择 嵌入 区 域 对 提高 水 印 抗 几何 攻击 能 力 至 关 重 要 。 原 因 


标识 ， 目 的 是 确认 内 容 创建 者 ， 判 断 载 体 是 否 被 自 改 。 水 印 算 ”是 水 印 图 像 受 几何 攻击 后 ， 载 体 图 像 和 水 印 坐 标 无 法 同步 ， 提 
法 的 核心 是 均衡 水 印 透明 性 与 鲁 棒 性 。 经 典 的 水 印 方法 有 奇异 ”取水 印 时 撕 入 区 域 不 能 准确 定位 ， 导 致 提取 水 印 不 完整 ， 因 此 
值 分 解 (SVD)、 正 交 三 角 分 解 (QR)、 三 角 分 解 (LU) 等 , 特别 是 奇 ”文献 [7] 首 次 提出 基于 特征 点 的 水 印 算法 。 文 献 [8] 利 用 SIFT 特 
异 值 分 解 有 效 提 高 了 水 印 算法 的 鲁 棒 性 。 文献 [1~3] 通 过 将 宿主 。 征 点 对 几何 攻击 不 变 特 性 ， 校 正 受 到 攻击 水 印 图 像 ， 恢 复 水 印 
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的 同步 性 ， 提 高 水 印 算法 抗 几何 攻击 能 力 ; 文献 [9] 将 SIFT 与 
DWT 结合 , 把 DWT 应 用 在 由 SIFT 选 定 的 特征 区 域 ; 文献 [10] 
利用 尺度 不 变 特 征 算法 提取 圆 形 区 域 特征 点 ， 将 水 印 典 入 在 
DCT 变换 后 的 中 频 系数 中 ， 抗 几何 攻击 能 力 较 强 ， 文 献 [11] 利 
用 SUREF 特征 点 构造 “和 鲁 棒 堆 水 印 ” 序 列 ， 提 高 水 印 算法 鲁 
性 。 文 献 [12，13] 将 SIFT 与 人 类 视觉 系统 (HVS) 相 结合 ， 根 据 
视觉 特性 确定 圣 入 强度 , 平衡 算法 的 透明 性 和 和 鲁 棒 性 。 文献 [14] 

将 果 蝇 优化 算法 (FOA) 应 用 在 小 波 域 的 水 印 算法 上 ， 自 适 
应 解决 了 透明 性 和 重 棒 性 的 均衡 问题 。 

综 上 所 述 ， 水 印 算 法 鲁 棒 性 与 透明 性 均衡 问题 一 直 在 发 展 
和 改进 中 , 但 仍然 存在 以 下 问题 a) 线性 尺度 空间 算法 不 能 准 
确定 位 特征 区 域 ;b) 线 性 构建 的 尺度 空间 提取 特征 点 不 够 稳定 ; 
c) 线性 尺度 空间 对 边缘 光滑 目标 提取 特征 点 不 精确 。 

本 文 提出 一 种 基于 非 线 性 尺度 空间 的 KAZE 算法 提取 特征 
点 ， 构 建 非 线 性 尺度 空间 不 变 特征 区 域 ， 以 增强 水 印 鲁 棒 性 ， 
通过 FOA 自 适 应 调整 泪 入 强度 ， 平 衡水 印 的 透明 性 和 重 棒 性 ， 
利用 非 线性 算法 选取 特征 点 定位 水 印 位 置 ， 再 做 水 印信 息 栓 入 
的 道 过 程 ， 提 取水 印 。 本 算法 相 比 于 线性 水 印 算法 ， 嵌 入 区 域 
定位 更 精确 ， 提 取 特 征 点 更 稳定 ， 更 能 有 效 抵抗 噪声 攻击 、 常 
见 攻击 和 几何 攻击 ， 自 适应 优化 算法 克服 了 非 线性 特征 区 域 构 
建 及 骨 入 强度 参数 随机 选取 工作 量 大 的 不 足 。 


1 ”基础 理论 


1.1 非 线 性 尺度 空间 构建 

非 线 性 尺度 空间 的 构建 主要 基于 非 线 性 扩散 滤波 原理 [5]， 
特点 是 通过 非 线性 扩散 ， 在 非 线 性 尺度 空间 检测 和 描述 特征 ， 
平滑 噪声 , 尽 可 能 保留 目标 原本 边界 ,使 其 特征 点 更 具 稳 定性 ， 
可 以 极 大 地 保持 图 像 的 细节 信息 。 非 线性 扩散 滤波 原理 与 物理 
学 中 二 维 热 扩散 过 程 描 述 相 似 ， 可 以 将 图 像 亮度 视 为 能 量 ， 亮 
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素 为 单位 的 尺度 参数 oi 转换 至 时 间 为 单位 如 式 (3〉 所 示 。 


50 G) 


KAZE 算法 构建 尺度 空间 的 流程 是 ， 首 先 通 过 高 斯 卷 积 平 
滑 处 理 图 像 ， 减 小 噪声 和 形变 等 干扰 ， 然 后 计算 图 像 梯 度 直方 
图 ， 以 获取 对 比 度 因子 参数 k， 确 定 保留 边缘 信息 量 ， 最 后 用 
进化 时 间 序 列 构造 非 线性 尺度 空间 ， 并 且 通 过 AOS 算法 及 简 
单 迭 代 求 解 ， 获 取 非 线性 尺度 空间 图 像 序列 。 如 式 〈4) 所 示 。 


1 
mm 


En -( 一 (bu -5) 立 4(O] 于 (4) 


I1=1 


其 中 : Al 是 i 尺度 上 图 像 亮度 Li 在 维度 1 上 的 传导 和 矩阵 ,I 为 单 
位 矩阵 。 
此 得 到 的 特征 点 有 更 高 的 重复 性 和 特殊 性 ， 能 提取 出 图 
像 中 鲁 棒 性 强 的 特征 点 ， 更 加 准确 地 确定 嵌入 区 域 。 
1.2 水 印 贬 入 与 提取 

SVD 相 比 于 其 他 水 印 嵌 入 方式 稳定 性 更 好 , 对 非 线性 特征 
区 域 做 奇异 值 分 解 ， 产生 两 个 奇异 向 量 矩 阵 (U、V) 和 一 个 对 
角 和 矩阵 (S) 。 奇 异 向 量 表征 图 像 结构 信 息 ， 对 角 线 元 素 值 包含 
图 像 重 要 特征 信息 。 图 像 受 到 一 定 程 度 干扰 时 ， 对 角 线 上 奇异 
值 
而 


ul 


不 会 有 太 大 波动 ， 仍 能 从 分 解 的 奇异 值 中 提取 出 水 印信 息 


量 由 于 奇异 值 是 矩阵 元 素 之 间 的 一 种 表示 03， 而 非 视觉 特性 
的 体现 ， 能 很 好 的 保证 水 印 的 不 可 见 性 ， 因 此 将 水 印 嵌 奇异 值 


根据 人 类 视觉 系统 特性 ， 选 择 适 当 的 嵌入 位 置 和 强度 可 以 
有 效 均衡 水 印 。 小 波 变换 多 分 辨 率 特性 和 人 眼 视觉 对 低频 分 量 
敏感 的 特性 一 致 ， 而 低频 分 量 又 聚集 大 量 能 量 信息 ， 因 此 将 水 
印 炙 加 在 小 波 分 解 后 的 低频 分 量 中 。DWT-SVD 结合 算法 使 水 
印 具有 更 好 的 透明 性 和 和 鲁 棒 性 。 
传统 水 印 幅 入 方法 是 通过 人 工 多 次 选取 幅 入 强度 参数 ， 工 


度 随 尺度 变化 模拟 热传导 流动 ， 通 过 热流 扩散 过 程 函 数 描述 图 
像 亮 度 扩散 。 该 函数 的 引入 使 得 扩散 能 够 适应 图 像 局 部 特征 ， 
更 准确 的 提取 出 图 像 特征 。 如 式 〈1) 所 示 。 


世 =div(c(x, y,1) VD) (1) 


其 中 : div 和 V 分 别 表 示 散 度 和 梯度 ， 函 数 c(x,y,t) 表 示 扩 散 的 
传导 函数 ， 参 数 t 是 尺度 参数 ，t 越 大 ， 图 像 表 示 形 式 越 简单 。 

KAZE 是 最 新 的 非 线 性 尺度 空间 算法 ， 它 通过 加 性 算 子 分 
裂 算法 (AOS) 和 可 变 传导 扩散 方法 来 构造 ， 并 进行 非 线 性 扩散 ， 
通过 任意 步 长 构建 稳定 的 非 线性 尺度 空间 09。 尺 度 级 别 按 对 数 
递增 ， 共 有 O 组 octaves， 每 个 octave 有 S 个 sub-level， 各 个 
层级 均 采 用 与 原始 图 像 相 同 的 分 辨 率 。 尺度 参 数 c 如 式 (2) 所 
不 。 


a.(0s)=0,2" ,0E[0...0-1],sE[0..S-1JiE[O0.N] ©) 


其 中 : ae 是 尺度 参数 的 初始 基准 值 , N=O*S 是 整个 尺度 空间 包 


作 量 大 旦 不 能 批 处 理 。 果 蝇 优 化 算法 能 通过 自 适应 迭代 ， 选 取 
水 印 嵌 入 的 最 佳 位 置 和 强度 ， 在 该 算法 中 ， 食 物 位 置 代表 嵌 入 
位 置 ， 味 道 浓度 值 代表 和 入 强度 。 该 算法 包括 以 下 几 个 步骤: 
a) 初始 化 最 大 迭代 次 数 Max 、 种 群 规模 Pop ， 随 机 初始 
化 果 蝇 群体 位 置 (x_fly,y_fly) 以 及 果 蝇 个 体 浓度 Sml 。 
b) 赋予 果 蝇 个 体 利 用 嗅觉 搜寻 食物 随机 方向 与 距离 ，Ram 
为 搜寻 距离 。 如 式 (5〉 所 示 。 


X_fly=x_fy(i)+Ram 
| 三 jp=y Ap(DHRam ( 5 ) 


c) 预计 与 原点 距离 为 D, 味 道 浓度 值 Si 为 距离 D 的 倒数 。 
如 式 〈6) 所 示 。 


D=sqrt(X_ fh (i) 2+Y fy(i) 2) 


S,; es (6) 


I~ 


d) 将 味道 浓度 值 S 代入 味道 浓度 适应 度 函数 ， 求 得 果 蝇 个 
味道 浓度 Sml(i)， 如 (7〉 所 示 。 


体位 


含 的 图 像 总 数 。 非 线性 扩散 滤波 模型 是 以 时 间 为 单位 的 ， 将 像 


Sml(i)=obj(S (i) (7) 
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适应 度 鞭 


数 影响 优化 速度 和 最 优 解 ， 即 兼顾 水 印 透明 性 和 


鲁 标 


其 中 : PSNR 是 所 


评价 
如 式 (8) 所 示 。 


性 。 采 用 峰值 信 噪 比 (PSNR) 和 归 一 化 系数 (NC) 作为 


标准 ， 分 配 PSNR 和 NC 的 权 值 比 ， 构 成 目标 函数 Obj。 


on=axlpsna -psn (1-o) -tSnc | 
na 


(8) 


期 望 的 峰值 信 噪 比 ，C 是 给 峰值 信 噪 比分 配 


的 权重 ，n 为 对 图 像 攻击 的 次 数 。 


e) 找 出 果 晶 群体 中 味道 浓度 最 高 的 位 置 《 最 优 值 ) ， 如 式 
(9) 所 示 。 
[bestSmell bestIndex| = min(Smel) (9) 
Smellbest = bestSmell 
保留 最 佳 味 道 浓度 值 和 与 其 对 应 的 位 置 坐标 (X_fly， 
Y_fly〉， 如 式 (10) 所 示 。 
X fly= X(bestIndex) (10) 
Y fly= Y(bestIndex) 
是 达 代 寻 优 ， 重复 执行 步 又 b)~e)， 并 判断 味道 浓度 是 否 优 


于 前 一 迭代 味道 浓度 ， 若 是 则 执行 步骤 用， 秋 代 次 数 达 到 最 大 
值 时 结束 循环 ， 得 到 最 佳 幅 入 参数 。 
2 ”算法 实现 
2.1 水 印 图 像 预 处 理 
为 消除 图 像 像素 空间 相关 性 , 需要 对 水 印 图 像 做 加 密 处 理 。 

采用 Anorld 对 水 印 图 像 中 各 像素 点 的 位 置 进行 随机 而 均匀 的 
置 乱 ， 从 而 消除 空间 像素 的 相关 性 。 置 乱 具 有 周期 性 ， 多 次 变 
换 图 像 会 在 某 一 时 刻 恢复 到 原来 的 图 像 ， 其 周期 与 图 像 的 像素 
大 小 有 关 。 图 1(a) 是 32X32 的 二 值 水 印 , 对 该 图 像 进 行 Anorld 
乱 ， 图 Cs 图 (d) 中 水 
像 在 和 24 时 ， 


2.2 载体 图 像 戏 入 位 置 


得 大 量 的 特征 点 ， 但 并 不 是 所 有 特征 点 都 适用 于 水 印信 息 


入 和 


a) 原 始 水 印 


特征 区 域 的 大 


提取 。 因 此 在 


除 掉 对 比 度 低 和 不 


定性 


在 相信 
过 计 
水 印 
取 对 


提取 处 理 ， 


。 筛 选 的 步骤 为 : 首先 采用 
法 ， 通 过 非 线性 滤 ; 


b) t=5 
图 1 不 同 次 数 的 水 印 置 乱 


c) {=20 d) t=24 


A 


水 印 图 像 大 小 确定 。KAZE 算法 虽 能 获 


嵌入 操作 之 前 ， 必 须 对 特征 
稳定 的 特征 点 ， 让 所 剩 特征 点 都 具有 较 高 稳 
KAZE 作为 图 像 特征 点 提取 方 
对 非 线性 尺度 空间 的 原始 图 像 进行 特征 点 


正点 进行 筛选 ， 删 


于 


闷 


邻 层 并 且 距 离 


算 生成 距离 集 


图 像 的 大 小 ， 


获得 特征 


点 的 集 {P} ,删除 掉 3x3 矩形 空间 内 两 点 
相等 的 重复 特征 点 ， 并 保存 为 {P") ， 然 后 通 
{D} ， 并 对 集 {D} 从 小 到 大 排序 ， 最 后 根据 
选择 特征 点 的 个 数 ， 按 照 集 {D} 中 的 顺序 选 


应 的 特征 点 和 


特征 区 域 作为 水 印 嵌 入 子 块 。 图 2 是 对 载体 


齐 向 明 ， 


图 像 Lena、Baboon、Barbara、plane 和 Camer 使 用 KAZE 算法 
点 较 均匀 地 分 布 在 载体 图 像 上 。 


提取 的 特征 点 ， 可 以 看 出 特征 


Baboon Barbara 


Plane Camer 
图 2 图 像 特征 点 的 提取 
2.3 自 适应 均衡 水 印 
2.3.1 均衡 水 印 说 入 


设 原始 载体 图 像 了 的 大 小 为 MxM， 水 印 图 像 到 的 大 小 
水 印 嵌 入 的 具体 流程 如 下 过 程 ; 
a) 根 据 2.2 节 提 取出 原始 载体 图 像 的 特征 点 ， 得 到 一 维 


构造 成 与 水 印 图 像 大 小 相同 的 矩阵 TT。 并 对 TT 人 

级 离散 小 波 变换 ， 取 三 级 低频 子 带 LL3 进行 SVD 分 解 得 : 
[U.S,V] = SVD(LL3) 

b) 将 置 乱 后 的 水 印 图 像 进 行 奇异 值 分 解 ， 得 到 Uw，Vw 

对 角 阵 Sw, 提取 右 奇 异 和 矩阵 和 对 角 和 矩阵 , 将 该 乘积 作为 水 印 

入 的 信息 ， 左 奇异 矩阵 用 于 水 印 的 提取 。 水 印信 息 和 矩阵 : 


为 


(11) 


和 
EA 


Waer=Sw (i )* Vy (i ) (12) 
oO) 使 用 式 (13) 在 S 中 添加 水 印 ， 得 骨 入 水 印 和 矩阵 WW， 式 

中 SSG) 由 果 晶 优化 算法 途 代 得 出 : 
WW =S +SS (ii)*W,., (13) 


dj) 对 WW 进行 奇异 值 分 解 ， 得 奇异 对 角 阵 S1， 将 其 与 


1) 


中 的 志和 严 使 用 式 (14) 得 峙 入 水 印 后 的 LL3 子 带 和 矩阵 CW。 
CW =U*S1*V” (14) 
6) 对 CW 做 小 波 逆 变换 ， 重 构 水 印 图像 。 对 含有 水 印信 息 

的 矩阵 降 维 ， 蔡 换 原始 图 像 的 特征 点 ， 得 到 水 印 能 入 图 像 。 
为 确保 算法 透明 性 ， 采 用 PSNR 作为 水 印 图 像 透 明 性 的 评 

价 标准 ，PSNR 值 的 定义 如 式 (15) 所 示 。 
eto Ga | (15) 
MSE 


其 中 : MSE 为 原 图 像 与 嵌入 水 印 后 图 像 的 均 方 误差 。 
2.3.2 水 印 提取 
a) 将 舱 入 水 印 图 像 的 相应 
使 其 构造 成 与 水 印 大 小 一 致 的 矩阵 ， 作 为 水 印 的 提取 位 置 。 
b) 将 该 矩阵 做 三 级 离散 小 波 变换 , 得 到 低频 子 带 LL3”。 
LL3” 做 奇异 值 分 解 ， 得 对 角 和 矩阵 S2 。 
c 将 b) 中 的 对 
和 V1 进行 奇异 


像素 点 提取 出 来 ， 呈 一 维 序列 ， 


对 


9 奇异 值 矩 阵 S2 和 内 入 水 印 时 的 和 珑 阵 Ul 
值 道 变换 ,结果 为 扎 阵 SN, 拘 入 强度 参数 SS(i) 


人 810.00077V1 


chinaXiv 


录用 定稿 
1.3 节 果 晶 优 化 算法 得 出 ， 使 用 式 (16) 提 取出 水 印 图 像 信息 
WN 矩阵 。 
WN=(SN-S)/SS(ii) (16) 
d) 利 用 水 印 的 左 奇异 矩阵 得 到 置 乱 的 水 印 。 
W’ =U, xWN (17) 
6) 使 用 Anorld 逆 变 换 得 到 原始 水 印 图 像 。 
采用 相似 度 评 价 方法 验证 算法 实用 性 ， 通 过 计算 原始 水 印 


WV 和 提取 出 水 印 丈 ' 的 相似 度 来 衡量 提取 水 印 的 有 效 性 ， 相 似 
NC 的 定义 公式 为 


2 WG DxW',D)) 


NS (18) 
> Gi, Dx Pw (i,)) 
其 中 : W(1)j) 为 原始 图 像 像 素 值 ，W “(i,j) 为 骨 入 水 印 后 图 像 的 像 


素 值 ，N 为 图 像 的 大 小 ，i j 为 矩阵 的 i 行 j 列 。 
3 “仿真 实验 与 分 析 


水 印 透明 性 用 峰值 信 噪 比 (PSNR) 衡量 , 值 越 大 水 印 透明 
性 越 好 ， 和 鲁 棒 性 用 归 一 化 相关 系数 (NO) 评价 ， 值 越 大 表示 和 鲁 棒 
性 越 强 ,选取 512X512 标准 灰 度 图 像 Lena、 Baboon、Barbara、 
Plane 和 Camer 作为 原始 载体 图 像 如 图 3(a)-(e), 选取 32X32 的 
二 值 图 像 “ 辽 宁 工 大 ”和 64 X64 的 二 值 图 像 “版权 保护 ”为 水 
印 图 如 图 3( 用 -(g)。 测 试 试验 中 KAZE 算法 的 Octave=2，Sub- 
level=4， 果 蝇 优 化 算法 种 群 规模 设 定 为 5， 最 大 和 迭代 次 数 为 50， 
PSNR 设 定 为 42， 对 水 印 大 小 为 32X32 的 “辽宁 工大 ”C 值 
设 定 为 0.4, 水 印 大 小 为 64X64 的 “版 权 保护 ”C 值 定 为 0.3。 
实验 中 分 别 对 水 印 载体 图 像 进行 以 下 9 种 攻击 : 压缩 图 像 ， 高 
斯 噪声 ， 乘 性 噪声 ， 椒 盐 噪 声 ， 中 值 滤波 ， 缩 放 图 像 ， 旋 转 图 
像 ， 剪 切 图 像 ， 平 移 含 水 印 图 像 。 


. 


I Et 


(a) (9) (d) 
辽宁 
工大 
保护 
(9 (@9) (g) 
图 3 原始 图 像 和 水 印 图 像 
3.1 透明 性 实验 
从 载体 图 像 提取 特征 点 时 存在 差异 ,因此 从 多 组 实验 中 (50 


次 ) 选取 一 组 最 优 解 。 图 4(a) 为 标准 灰 度 图 像 Lena、Baboon、 
Barbara Plane 和 Camer。 图 像 下 面 标 有 水 印 的 嵌入 强度 ss (ii) 
值 , 用 果 蝇 优化 算法 对 水 印 图 像 Lena、Baboon、Barbara、 Plane 
和 Camer 计算 最 优 仍 入 强度 ， 根 据 数据 可 以 分 析出 ， 使 用 不 同 
的 原始 图 像 和 水 印 图 像 所 得 到 的 最 优 和 入 参数 各 不 相同 。 图 4(b) 


为 嵌入 水 印 后 的 载体 图 像 ， 


值 。 


性 和 和 鲁 棒 


图 4(0) 为 提取 出 的 水 印 


ss(ii)= 1.39759 
本 Ta 


ss(ii) = 0.66449 


(a) 原 始 图 像 


Chi 


了 : 一 种 非 线性 尺度 空 


图 下 标 有 图 像 的 峰值 信 噪 比 PSNR 


区 


图 下 标 有 实验 得 到 的 NC 值 。 


PSNR=44. 1152 
(人 b) 含 水 印 图 像 


NC=1.0000 


NC=1.0000 


NC=1.0000 


NC=1.0000 


NC=1.0000 


版 权 


(oc) 水 印 图 像 


图 4 原始 图 像 .含水 印 图 像 和 提取 


8 的 水 印 图 像 


仿真 实验 表明 ， 本 文 算法 在 没有 受到 攻击 的 条 件 下 ， 透 明 


性 都 达到 了 理想 效果 。 从 结果 可 以 看 出 ，PSNR 值 都 
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达到 44dB 以 上 , 具 
时 ， 提 取出 来 的 水 印 NC 值 都 达到 了 1.0000， 
棒 性 。 


实验 结果 验证 
间 的 矛盾 。 
3.2 ”和 鲁 棒 性 实验 

为 检验 本 文 算法 抵抗 外 部 攻击 的 能 
Baboon、Barbara、Plane 和 Camer 五 幅 
击 和 几何 攻击 ， 验 证 其 鲁 棱 性 。 
3.2.1 常规 攻击 


了 较 好 的 感官 效果 , 水 印 的 透明 性 较 好 。 同 


水 印 具 有 较 强 外 


GhinaXiy 会 作 期 刊 ， 


齐 向 明 ， ee 性 尺度 空间 自 适 应 均衡 


规 攻击 后 含 64X64 水 印 图 像 和 提取 出 水 印 图 像 。 从 视觉 上 来 
看 ,所 提取 出 的 水 印 清晰 可 见 ,说 明 本 算法 可 以 抵抗 常规 攻击 ， 


从 视觉 效果 来 看 ， 含 水 印 图 像 与 原始 水 印 图 像 没 有 明显 
变化 , 图 像 质量 良好 , 且 提取 出 的 水 印 清晰 可 辩 ， 
了 本 文 算 法 能 较 好 的 均衡 水 印 透 明 性 


还 原 度 很 高 。 
与 鲁 棒 性 之 


， 分 别 对 Lena、 
含水 印 图 像 进行 常规 攻 


分 别 对 含 32X32 和 64X64 水 印 的 图 像 依次 进行 不 同 程度 


的 额 的 高 斯 噪声 、 乘 性 


噪声 、 椒 盐 噪 声 、 中 值 滤波 和 压缩 攻击 ， 


提取 水 印 NC 值 如 表 1、2 所 示 。 
表 1 含 32x32 水 印 图 像 常规 攻击 提取 NC 值 
攻击 方式 ”参数 ”Lena Baboon Barbara “Plane Camer 
We 0.01 0.9981 0.9985 0.9983 0.9979 0.9971 
高 斯 噪声 
0.02 0.9927 0.9932 0.9919 0.9916 0.9910 
0.01 0.9993 0.9989 0.9983 0.9992 0.9990 
乘 性 噪声 
0.02 0.9930 0.9928 0.9922 0.9931 0.9934 
0.01 0.9985 0.9968 0.9964 0.9956 0.9958 
椒盐 噪 
0.02 0.9921 0.9896 0.9911 0.9912 ”0.9904 
中 值 滤 波 3X3 0.9981 0.9926 0.9920 “0.9990 0.9966 
20 0.9999 0.9998 0.9999 1.0000 1.0000 
压缩 攻击 50 1.0000 1.0000 1.0000 1.0000 1.0000 
70 1.0000 1.0000 1.0000 1.0000 1.0000 
表 2 含 64x64 水 印 图 像 常规 攻击 提取 NC 值 
TT 攻击 方式 ”参数 Lena Baboon Barbara Plane Camer 
ll 0.01 0.9989 0.9982 0.9987 0.9989 0.9981 
高 斯 噪声 
0.02 0.9940 0.9958 0.9944 0.9943 0.9938 
0.01 0.9996 0.9991 0.9995 0.9996 0.9993 
乘 性 噪声 
0.02 0.9967 0.9957 0.9969 0.9963 0.9967 
0.01 0.9990 0.9971 0.9970 0.9975 0.9976 
椒盐 噪 
0.02 0.9923 0.9929 0.9903 0.9913 0.9909 
中 值 滤波 3X3 0.9990 0.9948 0.9935 0.9994 0.9983 
20 0.9999 1.0000 0.9999 0.9992 0.9996 
压缩 攻击 50 1.0000 1.0000 1.0000 0.9999 0.9999 
70 1.0000 1.0000 1.0000 1.0000 1.0000 
数据 表明 ， 随 着 攻击 强度 的 增 大 ， 对 提取 出 的 水 印 NC 值 
有 一 定 的 影响 ， 但 提取 出 的 水 印 NC 值 仍然 较 大 ， 特 别 是 受到 
噪声 攻击 和 压缩 攻击 ，NC 值 还 能 够 达到 0.99 以 上 ， 压 缩 程度 


到 达 70 时 , 所 有 图 像 NC 值 都 能 达到 1.0000， 
提取 的 质量 影响 不 大 。 对 比 表 1 和 表 2 


说 明 攻 击 对 水 印 


的 数据 , 嵌入 32X32 大 
小 的 水 印 和 64X64 大 小 水 印 在 遭受 相同 常规 攻击 时 得 


出 的 水 


印 效果 相似 ，NC 值 的 走向 大 致 相同 。 


图 5 


展示 了 常规 攻击 后 


含 32X32 水 印 图 像 5(a) 和 提取 出 水 印 图 像 5(b), 图 6 展示 了 常 


鲁 棒 性 较 强 。 


a(5)JPEG 压缩 (50)a(6) 高 斯 (10)a(7) 椒 盐 (0.02)a(8) 乘 性 (0.02) 

EAE 

工 E 工大 工大 工大 工大 I 

bD bb bb bb bg b(8) 
5 ”常规 攻击 图 像 和 32x32 水 印 图 像 


ee 


b(7) 


c(1) 高 斯 (0.01) c(O) 乘 性 〈0.02) ”cl3) 压 缩 (50) c(4) 压 缩 (70) 


版 权 畏 版 权 国 版 权 国 版 权 
保护 辆 保护 辆 保护 辆 保护 


d(1) d(2) d(3) d(4) 
图 6 常规 攻击 图 像 和 64x64 水 印 图 像 

3.2,2 几何 玻 击 

分 别 对 Lena、Baboon、Barbara、Plane 和 Camer 五 幅 含 水 
印 图 像 进行 不 同 程度 几何 攻击 。 表 3、 表 4 分 别 列 出 了 舱 入 32 
X32 和 64X64 水 印 图 像 的 NC 值 。 数 据 表 明 ， 旋 转 、 剪 切 和 
平移 攻击 后 ， 在 一 定 强度 内 ， 随 着 攻击 参数 强度 增 大 ， 对 图 像 
NC 值 产生 的 影响 不 大 。 嵌 入 32X32 和 64X64 大 小 水 印 图 像 
数据 的 走势 相同 , 而 且 柑 入 64X64 水 印 图 像 的 鲁 棒 性 更 优 , 这 
是 由 于 原始 图 像 的 特征 点 提取 数量 是 与 水 印 大 小 成 正比 的 ， 特 
征 点 提取 的 越 多 ， 在 受到 相同 攻击 时 ， 对 水 印 图 像 像素 值 的 影 
响 就 会 越 小 。 攻 击 实验 提取 出 的 水 印 NC 值 都 大 于 0.96， 说 明 
本 文 算法 可 以 较 好 地 抵抗 几何 攻击 。 


表 3 含 32x32 水 印 图 像 几何 攻击 提取 NC 值 
攻击 方式 参数 Lena Baboon Barbara Plane Camer 
、 10 0.9930 0.9828 0.9824 0.9839 0.9828 
- 30 0.9864 0.9818 0.9801 0.9676 0.9642 
剪 切 左上 1/4 0.9884 0.9809 0.9793 0.9751 0.9782 
ee (10,10) 0.9908 0.9848 0.9860 0.9862 0.9831 
(25,25) 0.9718 0.9705 0.9712 0.9710 0.9649 
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录用 定稿 齐 向 明 ， 等 : i NinaXiyE 全 信 明 , 
表 4 含 64x64 水 印 图 像 几何 攻击 提取 NC 值 表 5 本 文 算 法 与 文献 [10,14,18] 对 比 
攻击 方式 参数 Lena Baboon Barbara Plane Camer 攻击 方式 参数 文献 [10] ”文献 [14] ”文献 [18] ”本 文 算 法 
| 10 0.9987 0.9941 0.9940 0.9924 0.9928 ,0.01 0.9431 0.9819 0.9910 0.9981 
30 0.9885 0.9835 0.9837 0.9776 0.9772 I 00 0.9273 0.9774 0.9898 0.9927 
前 切 左上 1/4 0.9912 0.9819 0.9771 0.9766 0.9793 i 0.01 0.9429 0.9871 0.9958 0.9985 
M (10,10) 0.9916 0.9889 0.9884 0.9897 0.9867 0.02 0.9322 0.9801 0.9912 0.9921 
4 (25,25) 0.9872 0.9741 0.9737 0.9741 0.9720 乘 性 噪声 。 ”0.01 0.9329 0.9814 0.9949 0.9993 
从 下 面 展 示 的 部 分 含水 印 载体 图 像 和 提取 出 的 水 印 图 像 来 50 0.9317 0.9962 0.9887 1.0000 
看 ， 含 水 印 载体 图 像 质量 良好 ， 所 提取 出 的 水 印 与 原始 水 印 图 70 0.9392 0.9976 0.9909 1.0000 
像 并 无 差别 ， 水 印 图 像 清晰 可 见 ， 还 原 度 很 高 。 实 验 验 证 了 本 中 值 滤波 3X3 0.9740 0.9912 0.9578 0.9981 
文 算法 在 透明 性 良好 的 前 期 下 ， 水 印 仍 具有 较 强 鲁 棒 性 ， 达 到 | 10 0.8838 0.9704 0.9978 0.9930 
了 水 印 算 法 透明 性 和 和 鲁 棒 性 的 均衡 。 部 分 受 几 何 攻 击 的 图 像 如 i 30 0.8046 0.9320 0.9886 0.9862 
图 7、8 所 示 。 剪 切 左上 1/4 0.9504 0.9785 0.9875 0.9884 


平移 (10,10) 0.8692 0.9725 0.9425 0.9908 


xl x2 x3 x4 x5 x6 x7 x8 x9 xl0 xl xl2 
类 


攻击 类 型 


a(5) 旋 转 (30) a(6) 剪 切 左上 角 1/4 a(7) 平 移 (10,10〉a(8) 平 移 (25,25) 图 9” 归 一 化 相关 系数 对 比 
a 为 了 验证 本 文 算法 在 具有 较 强 鲁 棒 性 的 同时 能 兼顾 水 印 的 
工人 了 工头 工大 工大 工大 透明 性 ， 选 择 含 32X32 大 小 水 印 的 Lena 图 像 的 PSNR 值 与 文 
bl) bo bj bdh bs boO bo bg) 献 [10，14，18] 对 比 。 本 文 算法 、 文 献 [10，14，18] 的 PSNR 值 


图 7 几何 攻击 图 像 和 32x32 水 印 图 像 分 别 为 48.6883、35.79、43.2748 和 48.4723， 且 本 文 算 得 到 的 

二 所 有 含水 印 图 像 的 PSNR 值 都 在 44dB 以 上 , 优 于 文献 [10,14， 
18]， 显 示 出 本 文 算 法 得 到 的 水 印 图 像 透明 性 良好 。 本 文 算法 对 
水 印 答 入 区 域 定位 更 准确 ， 提 取出 的 特征 点 更 稳定 ， 对 各 种 攻 
击 抵抗 力 较 强 ,说明 本 文 算法 可 以 均衡 水 印 的 透明 性 和 和 鲁 棒 性 。 
c( 有 ) 旋 转 (10) ” ”c(@2) 旋 转 (10) c(3) 剪 切 左 上 角 1/4 c(4) 平 移 (25,25) 结束 语 
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d(1) d(2) d(3) d(4) 的 水 印 嵌入 区 域 选 取 方 案 ， 依 据 果 晶 砚 食 行为 寻求 自 适 应 均衡 

图 8 几何 攻击 图 像 和 64x64 水 印 图 像 水 印 最 优 参数 ， 结 合 离散 小 波 变换 和 奇异 值 分 解 算法 ， 增 强 水 

3.3 对比 实验 印 图 像 稳定 性 。 实 验 结果 表明 ， 将 水 印信 息 嵌 入 到 重要 特征 区 
以 Lena 为 原始 图 像 ， 32X32“ 辽 宁 工 大 ”为 水 印 图 像 , 实 。 域 中 ， 水 印 算法 鲁 棒 性 显著 提高 。 在 噪声 、 常 规 和 几何 攻击 条 


仿 结 果 与 文献 [10, 14, 18] 算 法 对 比 , 在 抵抗 噪声 和 常规 攻击 两 件 下 ， 与 现 有 文献 比较 ， 本 文 算法 在 水 印 均 衡 及 算法 效率 两 个 

个 方面 ,NC 值 都 达到 0.99 以 上 ,普遍 高 于 文献 [10，14，18]， 方面 都 有 提高 。 在 计算 速度 和 抗 旋转 攻击 方面 本 文 算法 还 存在 

说 明 本 文 算法 具有 较 强 鲁 棒 性 。 但 旋转 攻击 较 文 献 [18] 略 低 , 这 ” 不足， 有 待 进一步 研究 。 
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